hep-ph/9901323) 



IRB-TH-8/98 



Inclusive decays and lifetimes of doubly charmed 

baryons 

B. Guberinat B. Melicj H. Stefancic* 

Theoretical Physics Division, Rudjer Boskovic Institute, 
P.O.B. 1016, HR-10001 Zagreb, Croatia 



Abstract 

The analysis of singly charmed hadrons has been extended to the case of doubly charmed 
baryons, S+ and f2+. Doubly charmed baryons are described as a system containing 

a heavy cc-diquark and a light quark, similarly as in a heavy-light meson. This leads to 
preasymptotic effects in semileptonic and nonleptonic decays which are essentially propor- 
tional to the meson wave function. Interplay between preasymptotic effects in semileptonic 
and/or nonleptonic decay rates leads to very clear predictions for semileptonic branching 
ratios and lifetimes of doubly charmed baryons. 

PACS numbers: 12.38.Lg, 13.30.Ce, 13.30.Eg, 14.20.Lq 



1 Introduction 



Weak decays of heavy hadrons ||, E|, || H present a very rich field of phenomena owing to the 
complexity of confinement. Being the bound states of a heavy-quark and light-constituents 
(mesons, singly charmed/bottom baryons) or even of two heavy quarks and one light con- 
stituent (doubly charmed/bottom baryons), heavy hadrons contain soft degrees of freedom 
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which generate nonperturbative power corrections, such as the destructive and/or construc- 
tive Pauli interference and the W-exchange/annihilation between the quark coming from the 
heavy-quark decay and one of the light constituents. 

Inclusive decay rates and lifetimes of charmed mesons were reasonably reliably calculated in 
the last decade. The overall picture emerging is qualitatively satisfactory and the lifetime 
hierarchy has been predicted for singly charmed baryons and found to be in agreement with 
present experiments @. Also, the difference in lifetimes (a factor of 2-3) between D + and D° 
mesons, which is due to the negative Pauli interference preasymptotic effect, was explained a 
long time ago [|, |, [|. 

The numerical calculations performed in the middle of the eighties j|, |9| provided us with the 
predictions of a lifetime pattern which has recently been confirmed by experiment. This success 
is rather surprising since with the advent of the systematic Operator Product Expansion (OPE) 
[]l]] and Heavy Quark Effective Theory (HQET) [|l(]] it has become clear that the charmed quark 
mass is not heavy enough for the m" 1 expansion to be trustworthy. Nevertheless, it seems that 
employing systematically field-theory methods to the very end of the calculation (up to the 
hadronic wave function for which we have to rely upon some phenomenological models), one 
is able to make clear predictions to be compared with present and future experiments and 
possibly to disentangle between various preasymptotic effects. 

On the other hand, the inverse bottom quark mass appears to be a good expansion parameter 
in bottom decays. However, the role of four quark operators is negligible there (effects of the 
0(1%)) leaving charmed hadron decays as a playground for studying such effects and as a test 
of the possible violation of the quark-hadron duality. 



In this paper we extend the analysis of singly charmed baryons decays and lifetimes [11] to 
the case of doubly charmed baryons. Recently |fE^| , a rather phenomenological approach using 
effective constituent quark masses and fit of singly charmed baryon decays has been employed 
to study doubly charmed baryon decays. We , however, have used a systematic field-theory 
approach to the very end in order to be consistent with the previous treatment of singly charmed 
hadrons and have also included the preasymptotic effects in semileptonic decay rates of doubly 
charmed baryons, thus calculating all decay rates at the Cabibbo subleading level. We show 
that preasymptotic effects dramatically change the simple spectator picture and lead to a very 
clear pattern of semileptonic branching ratios and lifetimes. 



2 Preasymptotic effects and wave function in doubly charmed 
baryon decays 

Using the optical theorem, the inclusive decay width of a hadron H cc with mass Mh cc containing 
two heavy c quarks can be written as 

T(Hcc - /) = 777-^2 Im(H cc \f\H cc ) , (1) 



where T is the transition operator 
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f = ij d 4 xT{L eff (x),Ll ff (0)}. 



(2) 



In the following we use the Operator Product Expansion which is based on the assumption 
that the energy release in the decay of a c quark is large enough. This implies that momenta 
flowing through internal lines are also large and therefore justify the OPE. 

The general formula for the decay is given by (l|, ||. [| 



r(ir„-/) = ^\ y \ 2 ^}4{H cc \cc\H cc ) (3) 

f \Hcc\ 

+ C 5 9 

+ V 4 {H ^ Tiq) f ic) \ Hcc) +Q (l/mj)}. 
t— 1 mf. 

f f 

Here c s and c 5 are Wilson coefficient functions which are known at one-loop order and tree 
level, respectively Jl], |2|, ||]. V represents appropriate matrix elements of the CKM matrix. 

Let us calculate the semileptonic decay rates first. The main contribution is expected to come 
from the quark decay- type diagram, which is proportional to (H cc \cc\H cc ) and is given, to 
0(m 2 ), as 



rg(ff«) = 2-^mHch, SL (x)P (x) + s 2 „ SL (0)) 

v 2 m 2 2 m 2 7 v ' 

Here /j^(H cc ) and ijlq{H cc ) parametrize the matrix elements of the kinetic energy and the 
chromomagnetic operators, respectively. Their determination will be discussed later. 

Throughout the paper we use the abbreviations s 2 and c 2 for sin 2 9 c and cos 2 6 c (0 C is the 
Cabibbo angle). 

The next contribution is coming from the dimension-five operator 



r&(tf«) = 2-^3^(c 2 P l( x) + , 2 )(-2^££l ) . (5 ) 

Note that in both (g|) and (|) there is an additional factor 2 coming from the decays of two c 
quarks in the doubly charmed baryon. 

The phase-space corrections Pq and Pi are cited explicitly in the Appendix. The radiative 
QCD correction rjsi 0, |14| ^ s given by 



i \ 1 2a 5 

3 7T 



(6) 



where for g(x) we have 

25 

=7r 2 - — + x (18 + 87T 2 + 241nx), (7) 
and x = m 2 /m 2 . Leptons are taken to be massless. 

Recently, Voloshin has noticed [15] that preasymptotic effects in semileptonic inclusive decays 



can be very large owing to the constructive Pauli interference, the result up to the CKM matrix 
element being given by 



fsL = ^m c 2 (4v^-l)5|V>(0)| 2 . (8) 

Here k is a correction due to the hybrid renormalization of the effective Lagrangian and it takes 
care of the evolution of L e ft from m c down to the typical hadronic scale fi ~ 0.5 — 1 GeV . The 
factor 5 in front of |V>(0)| 2 reflects the spin structure of doubly charmed baryons. The baryon 
wave function vp(0) will be discussed later. 

The total semileptonic rate for one lepton species is given by 



T S l(H cc ) = Tf£(H cc ) + T G SL (H CC ) + T v s t° shm (H cc ) , (9) 

where 



-pi/ oloshm /"H — h\ q 

SL \ 'cc ) ' 

rV oloshin i c+ \ „2-p 
SL {-cc) = S L S L, 

rVoloshin^+j = c 2p 5L> 

(10) 



In view of the significant preasymptotic effects in the SL decay rates of singly charmed baryons, 
one can expect a large Pauli-interference contribution in the semileptonic decay rate of the f2+ 
baryon (cos quark structure), where that contribution is present at the leading Cabibbo level. 

Nonleptonic decay rates are slightly more complicated, since in the final state the lepton pair is 
substituted by a quark pair. The contributions analogous to (Q) and (|5|) are ( including 0(m^.) 
corrections) 



■pdec 
1 NL 



1927T 3 

„2„2\ 



m 5 c (c 2 _ + 2c 2 + )[({c 4 + s 4 )P (x) + 



c z s z )ri NL (x) + c 2 s 2 P (x)fi NL (x)] 

l /ll(Hcc) | l ^Hcc) 

2 m 2 2 m 2 . ' 



[1 



(11) 
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T% L (H CC ) = 2^m 5 c {(2cl+c 2 _)[((c 4 + s 4 )P 1 (x) + 

c 2 s 2 )t] NL (x) + c 2 s 2 P 1 (x)fj NL (x)\ 
+ 2(4 - c 2 _)[((c 4 + s A )P 2 {x) + c 2 s 2 ) VNL (x) + 

J2„2 



'•VA(. ( ) ; )vtir)!]}(~2^g^). (12) 



Radiative corrections to the nonleptonic decay, tjnl(x) and t\nl{x), are far more complicated 
than analogous corrections (g) and (^) to the semileptonic decay and the reader is referred to 
the original paper where they were first calculated |l6t| . 

Again, the preasymptotic effects are expected to contribute significantly to the total nonleptonic 
decay rate. They are given by 



r« = |L 2 [c 2 + |(i-^)(4-c 2 )]5|^(o)| 2 , 

= %rn 2 c [-\c + (2c„ - c + ) 

-~(1 - v^)(54 + c - " 6c + c -)] 5 I^(0)| 2 , 
r +' = f^[^ + (2c-+c + ) 

-hi - v^)(54 + + 6c+c_)] 5|^(0)| 2 . (13) 

An explicit calculation leads to the following nonleptonic decay rates: 



TNL{^tc~) — r^£(s+ + ) + r 7VL (r,+ + ) 

+ {(c 4 + s 4 )Pm*(z) + cV(l + P mt (x))}T M , 

("cc) = r^£(E+) + r^ L (H+) 
+(c 4 p e:c (x) + c 2 s 2 )r ea; 
+( s 4 p^(x) + c V)rf, 
r^(n+) = r#£(n+) + r£ L (n+) 
+( c 4 + cVp int (*))rf 

+(c 2 s 2 P ex (x) + s 4 )T ex . (14) 



All corrections P and P are given explicitly in the Appendix. 

An important remark to be made here concerns the mass parameters in the calculation of the 
matrix elements /j, 2 and \Xq. Whenever we perform an expansion which is essentially a field- 
theoretic procedure (either the OPE for the transition operator T or the HQET expansion in 
the case of the cc operator), the expansion parameter is always the current heavy-quark running 
mass m c . On the other hand, in the calculation of the matrix elements, which is performed 
within quark models, it is more appropriate to use constituent quark masses m* . 
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Following this procedure, we give the expressions for /i 2 and p? G . For fj, 2 , we have 



2 22 



"2m* 2 + m*m* 



+ 



2m] 



■)m 



2 

c ' 



(15) 



where u c is the average heavy-quark velocity in the ccq baryon, m* and m* are constituent 
masses of the heavy and the light quark, respectively, and T is the average kinetic energy of 
the light quark and the heavy diquark. The precise description of this calculation is given in 
[12] and relies upon some phenomenological features of the meson potential. 

The contributions to the /j,q operator, connected to the matrix element of the chromomag- 
netic operator, can be divided into two parts. The first part includes effects coming from the 
heavy-light chromomagnetic interaction and these contributions can also be found in the singly 
charmed baryon The second part comprises effects originating within the heavy diquark, 
i.e. heavy-heavy chromomagnetic interactions. These effects are new [12, 17] and characteris- 
tic of doubly charmed baryons. Their estimation relies upon the nonrelativistic QCD model 
calculation fl2], 17, 15, 19]. The final expression is 



= \{M* ccq - M ccq )m c - (§fl&^|L + \&^-) , (16) 

where the first term describes the heavy diquark-light quark hyperfine interaction, while the 
second and the third correspond to the interaction of two heavy c-quarks in a diquark state. 
They are of the "chromomagnetic" and "Darwin" type, respectively. In (|l6|), M ccq is the mass 
of the doubly charmed baryon, M* cq is the mass of its 3/2 spin counterpart and 0(0) is the 
wavefunction of the cc pair in the heavy diquark, i.e. \<p(0)\ 2 is the probability for these two 
heavy quarks to meet at one point. 

In the calculations above, up to the hadronic matrix elements we have used field theory only. 
The results are expressed in terms of the baryon wave function ^(0) and the matrix elements of 
the kinetic and chromomagnetic operators, which are /i 2 and Hq, respectively. The use of the 



usual singly charmed baryon wave function ^(O), as given in [11], would be premature, since 
intuitively, one expects a two-heavy-quark system to behave differently from the single-heavy- 
quark one. 

In the case of singly charmed baryons, the heavy quark is "sitting" in the centre of the baryon 
and the other two light quarks are moving around. Their spin and the color charges are 
correlated (in order to have the appropriate spin and color structure of the entire baryon) , but 
their spatial motion is not. In this way, one has a three-body picture of the baryon containing 
a single heavy quark and one should use the baryonic wave function ^(0) accordingly. 

In the case of doubly charmed baryons, one assumes that two heavy quarks are strongly bound 
into a color antitriplet state. As far as the light quark is concerned, the bound state of two 
cc-quarks appears as a pointlike diquark object |2^, 21, ^]. Thus, in the heavy-quark limit, 



which can (presumably) still be applied in our case (m c > Aqcd), a doubly charmed baryon 
appears to consist of a heavy diquark and a light quark forming a "meson" state. Therefore, 
one expects that the wave function of the doubly charmed baryon (which has to be considered 
as the light-quark wave function at the origin of the cc-diquark) behaves essentially as the 
mesonic wave function. 
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Figure 1: Dependence of nonleptonic decay widths on the wave function squared. The shaded 
area represents the unphysical region of negative Ttvl to stress the problem of choice of the wave 
function. The vertical line represents the value of the wave function squared (corresponding to 
/d = 170 MeV) used in the calculations. 



We use the derivation of hyper fine splittings of mesons calculated in the constituent nonrela- 



tivistic quark model by De Rujula et al |23|, 24] to obtain the following relation between the 



wave functions of the doubly charmed baryon and the D-meson: 

m ^ = l ]m?D = l!hMEf!l. (17) 

The factor 2/3 comes from the different spin content of doubly charmed baryons, i.e. the cc- 
diquark forms the spin-1 color antitriplet state. The baryonic wave function squared in ( |l7|) 
is directly proportional to the D-meson decay constant, fp, squared. The factor k -4 / 9 is the 
effect of the hybrid renormalization which accounts for the fact that fo is measured at the 
scale ~ m c (k = 1), and one has to evolve fo down to the hadronic scale /i = 0.5 — 1 GeV. 

The choice of the mesonic wave function ~ fz>, instead of the singly charmed baryon wave 
function ^(O)] 2 ~ Fp, where Fjj is the static value of the D-meson decay constant, also 
seems to be consistent numerically. In Figure 1. we have displayed the dependence of the 
Tnl(H cc ) on |?/>(0)| 2 in the large range of the fn values. In our numerical calculation we use 
fo = 170 MeV as a central value. This value is consistent both with QCD lattice calculations 



5fl and QCD sum rules calculations [26, P7[| . In the case of TnlC^Jc) there is a negative 



Pauli interference which cancels the contribution coming from the decay- type diagram (11) 



and the chromomagnetic operator (12). This case is very similar to the role of the negative 
Pauli interference in the D + decay where it competes with the decay-diagram contribution. For 
fo large enough, the nonleptonic and total rates in both D + and S+ + decays become negative, 
Fig.l. However, a reasonable choice of fo gives positive results. 

In view of these facts and results, we may conclude that the phenomenological rule of using 
fo in mesonic and its static value Fjj in baryonic systems, employed first in singly charmed 
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Figure 2: Dependence of lifetimes on the parameters \i and Aqcd- Both pictures show the 
evident insensitivity of T(fi+) and r(S+) under variation of fi and Aqcd, and a small but 
notable sensitivity of r(H+ + ). 



hadrons 0], can be successfully extended to the consideration of doubly charmed baryons. Not 
doing so, but taking Fjj instead of fn would lead us to the unphysical region. 

So, our result is a confirmation of the above mentioned phenomenological rule at the same, 
qualitative level. Taking this rule as postulated for singly charmed hadrons, we can interpret 
our results as an extension of the same rule into the doubly charmed sector. Also, we can 
generalize the rule to some extent. We see that the use of fo is required in singly charmed 



mesons and doubly charmed baryons, while Fjj is used in singly charmed baryons [11]. So, it is 
allowed to say that fo should be used in systems with two-body dynamics (in the ccq baryon 
case, a heavy diquark and a light quark) and Fd in systems with three-body dynamics. It is 
important to stress that these considerations and conclusions are of purely phenomenological 
origin, i.e. they have no direct justification in field theory. 



3 Semileptonic inclusive rates and lifetimes 
cussions 



results and dis- 



In numerical calculations we use the following set of parameters, which closely follows the set 
used in [11|. For Aqcd = 300 MeV, the Wilson coefficients are c + = 0.73 and c_ = 1.88. 
The charmed quark mass is taken to be m c = 1.35 GeV and for the strange quark mass we 
use m s = 150 MeV. The value of the average kinetic energy T, appearing explicitly in (|i~5|), 
is taken from [12] to be T = 0.4 GeV, and the light and heavy-quark constituent masses are 



in 



0.3 GeV and m* = 1.6 GeV, respectively. The numerical value of the diquark wave 
function is also taken from [12] to be |<^(0)| = 0.17 GeV 3 ^ 2 . The numerical values for masses of 
doubly charmed baryons are taken from 12T 



As far as the Aqcd dependence is concerned, in the range Aqcd ~ 200 — 300 MeV the lifetimes 
of H+. and 0+ are practically constant and the lifetime of is more sensitive to the value 
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^++ 






Nonleptonic widths in ps 1 




0.345 


4.158 


2.859 


Semileptonic widths in ps 1 




0.151 


0.173 


0.603 




Semileptonic branching ratios in 


% 


BRsl 


23.4 


3.9 


14.9 


Lifetimes in ps 


T 


1.55 


0.22 


0.25 



Table 1: Predictions for nonleptonic widths, semileptonic widths, semileptonic branching ratios 
(for one lepton species) and lifetimes of doubly charmed baryons for the values of parameters 
m c = 1.35GeV, n = lGeV, A QC d = 300MeV, f D = 170MeV. 
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Figure 3: Dependence of semileptonic decay widths and lifetimes on the value of the wave func- 
tion squared. The vertical line represents the | Y?(0) | 2 used in calculations which corresponds 
to fo = 170MeV. The second picture shows the instability of r(H+ + ) for large ^(O)] 2 . 



of Aqcd, being somewhat (10%) larger for Aqcd = 200 MeV. The same is true for the /Li- 
dependence in the reasonable range \i ~ 0.5 — 1 GeV. The lifetimes of E+ and 0+ stay almost 
constant with variation of /i, and the lifetime of H+ + grows slowly with fj, (by 18%), Fig. 2. 

From Table 1 one can see that the Voloshin type of preasymptotic corrections in the semilep- 
tonic decay rates of 0+ is significant, contributing at the Cabibbo leading level, Eq. (|To|). This 
contribution brings Tsi^tc) to be four times larger than TsL(^tc~)i which receives contri- 
butions only from @ and @ . In the TsiC^tc) there is the Pauli interference effect at the 
Cabibbo suppressed level, but it still makes the rate larger by 15% than that for the H+ + 
bar yon. 

Clearly, since both semileptonic and nonleptonic rates are significantly affected by the large 
preasymptotic effects which are proportional to |?/>(0)| 2 ~ /J,, the results for lifetimes and the 
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semileptonic branching ratio for fi+ depend crucially on the choice of fjj. The latter is obvious 
from Fig. 3, where r,sx(f2+) grows linearly with and especially from the second picture in 
Fig. 3, where it is clear that r(H+ + ) shows instability for fjj larger than 180 MeV owing to 
the large cancelation between the negative Pauli interference term and the contributions from 



Eqs. (11) and (|l2|). This is a clear signal that one should not take the results for too 



literally. 

Keeping the above remarks in mind, we predict the following pattern for semileptonic branching 
ratios: 



BR SL (E+) « BR SL (Sl+) « BR SL (E++) , (18) 
and the following pattern for the lifetimes: 

r(S+)~T(fi+)«r(H++). (19) 

We can compare our results with the calculations of lifetimes of H+ + and EjJ,, which appeared 
recently |l2|]. The authors of that paper employed a similar field-theory technique, but had 
a different approach to the choice of relevant parameters. Throughout their paper they used 
the constituent heavy-quark mass as an expansion parameter, which is a phenomenological 
procedure that we do not find fully consistent. In the calculation of semileptonic decay rates, 
they did not include large preasymptotic effects which significantly change total semileptonic 
widths. Comparison shows that our numerical results are significantly different from those of 
Kiselev et. al. ll! 



Finally, it is worth discussing briefly the decay of the heaviest weakly decaying charmed hadron, 
the triply charmed baryon Although its complicated structure and its intrinsic tree-body 

motion prevent us from applying to this particle the heavy-light picture as described above 
to other weakly decaying heavy hadrons, it is possible to give some qualitative predictions 
for the ^2cc d eca y ra te an d the lifetime. In this baryon, preasymptotic effects (giving large 
contributions in the singly and the doubly case) do not exist for lack of light valence quarks. 
Thus, the dominant contribution comes from the operators of dimensions three and five. Since 
in doubly charmed decays the contribution of dimension-five operators represents less than 20% 
of the contribution of the decay (dimension three) operator, it seems reasonable to approximate 
the total decay width of ^2++ with the triple c-quark decay contribution and estimate the error 
of disregarding dimension- five operators at the level of 20%. In this case, the expression for 
TTOT^tcc) can be obtained by multiplying the expressions @, (||), (|ll|) and (12) by a factor 



of 3/2, summing them and taking the limit ^ — > and Hq — ► 0. The numerical value for the 
lifetime is 



T(n+t) = 0A3ps. (20) 

As the calculation of dimension-five operator contributions in triply charmed baryon decay 
rates is out of scope of the present paper, this result can be considered only as qualitatively 
correct. 
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Figure 4: Dependence of lifetimes of weakly decaying charmed hadrons on their mass. Ex- 
perimental results are drawn as points with error bars (for singly charmed hadrons), while 
theoretical predictions are drawn as filled diamonds (all weakly decaying charmed baryons). 
Values for masses of doubly and triply charmed baryons are taken from p£| . 



4 Conclusions 



Application of the heavy-quark expansion to the problem of inclusive decays of doubly charmed 
baryons enables us to give very interesting predictions for their lifetimes and semileptonic 
branching ratios. Large lifetime differences are present between H+ + on the one hand and 
and fi+ on the other. Our numerical results pick out H+ + as the longest living charmed 
particle (Fig. 4), although the numerical value for r(H+ + ) should be taken with certain reserve 
for reasons already mentioned. Such a large numerical difference within the lifetime hierarchy 
makes these predictions suitable for testing by the first forthcoming experimental observation 
of doubly charmed baryons. A theoretical prediction for semileptonic branching ratios is even 
clearer and the hierarchy of BRsl is unambiguously determined. 

The total hierarchy of lifetimes for charmed hadrons is shown in Figure 4. It is evident that 
charmed hadrons show a very complex pattern in the r — M plane. One can note that doubly 
charmed baryon lifetimes are comparable with those of their singly charmed counterparts. This 
result is just opposite to the naive expectation of roughly double widths in the doubly charmed 
case owing to the decay of two, instead of one c quark, or, correspondingly twice smaller lifetimes 
of ccq baryons. However, the mesonic nature of doubly charmed baryon wave functions, where 
one uses the smaller fjj constant instead of its static Fd value, reduces four-quark operator 
contributions, increasing doubly charmed baryon lifetimes. 

Although the c quark is considered not being heavy enough to ensure a reasonable convergence 
of the heavy-quark expansion series, the numerical results in the singly charmed sector show 
satisfactory qualitative and even quantitative agreement with experiment [11]. If future exper- 



iments concerning the doubly (and triply) charmed sector should show similar agreement with 
our theoretical predictions, that would have important implications on the role of four-quark 
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operators and the entire theory of heavy quark expansion. Besides, the absence or the presence 
of agreement might have notable implications on the validity of the quark-hadron duality in 
the charmed sector. 



A Appendix: Phase-space corrections 

Phase-space corrections can be understood as reduction of particle phase space due to the 
propagation of the massive particle in the loops of diagrams describing inclusive decays. In our 
case, we consider the s quark to be massive, while the other particles (u,d quarks and leptons) 
are treated as massless. These corrections can be classified according to the type of the the 
operator diagram in which they appear and according to the number of massive quarks in the 
loop (or in a final state if we consider an inclusive process as a sum of exclusive channels) . 

First, we shall enumerate corrections that appear in decay and dimension-five operator dia- 
grams. From here on, x = m 2 /m 2 . 

• One massive quark in the loop: 

P (x) = (1 -x 2 ){\ - %x + x 2 ) - 12x 2 In x, (21) 

P l {x) = {l-x)\ (22) 

P 2 (x) = (1-xf. (23) 

Pq{x) appears as a correction to the decay- type diagram, while P\(x) and P2(x) come as 
corrections to the chromomagnetic operator. 

• Two massive quarks in the loop: 
Using the notation 



we have 



,{x) = y/l-Ax, (24) 



Po(x) = v(x)(l - Ux - 2x 2 - 12x 3 ) + 

24x 2 (l-x 2 )ln( 1 -±^\), (25) 
1 — v(x) 



Pi(x) = -(2P (x) - ydyPo(y) \ y=x ) , 



P2(x) 



■> x . , 9n , , \ + v(x) . 
v(x)(3x 2 + - + 1) - 3x(l - 2x 2 ) In ( 



1 — v(x) ' 



(26) 
(27) 



12 



Similarly as above, Pq(x) appears in the decay diagram, while P\(x) and P2(x) are corrections 
to the dimension-five operator. 

Corrections in the paper due to one massive quark are systematically denoted by P, while those 
due to two massive quarks are denoted by P. 

Next, we display the phase-space corrections to four-quark operators: 

P ex (x) = (1-x) 2 , (28) 

P int (x) = (l-x) 2 (l + x), (29) 



P int (x) = VT^4x~. (30) 

P ex (x) appears as a correction to the exchange diagram, Pi n t(x) corrects for the massive quark 
in the interference contributions, while Pi n t{x) is a correction in the case of negative interference 



when there are two massive quarks in the loop y, 3C, 31 1. 
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